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, Abstract. We report on recent results obtained for the massive operator matrix elements which 

£H ■ contribute to the massive Wilson coefficients in deep-inelastic scattering for Q 2 ^> m 2 in case of 
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S : INTRODUCTION 

^! To determine the value of the strong coupling constant a s (M|) from the world deep- 

inelastic scattering data [1] precisely, the structure function F2(x, Q 2 ) has to be described 
to O(cCs) accuracy at present. The corresponding next-to-next-to-leading order (NNLO) 
t-h ■ parton distributions are instrumental for accurate predictions of the weak boson and 

^ ! Higgs boson production cross sections at hadron colliders J20. While the massless 

\ contributions are known to NNLO yfl, in the massive case this has been achieved for 

a series of Mellin moments N only in the region Q 2 > lOm^ 0] 0, with mg the mass 
of the heavy quark. Based on the 2-loop massive operator matrix elements, which are 
known for general values of N [Q] to the linear term in the dimensional parameter e 
||s|], and the massless 3-loop Wilson coefficients [@] one may compute all logarithmic 
contributions at 3-loop order iTToh . In the kinematic region of HERA, however, these 
j> | terms alone are not dominant over the yet unknown constant term for general values of 

N, as has been demonstrated for the moments calculated in ||4j]. The O(a^) corrections 
to the structure function Fl(x, Q 2 ) in the asymptotic region were calculated in (LI], 3]. 
Heavy flavor corrections for charged current reactions were given in IU2I1 . 

The heavy flavor contributions to the structure function F2(x, Q 2 ) at O(a^) in case of 
one massive quark are described by the five massive Wilson coefficients in the asymp- 
totic region L^ s ,Z^ S ,lj,#£ s 0]. Two of these Wilson coefficients, L^ s and Zj, 



have been computed completely for general values of Af in 111311 . cf. also 111 Oil . In II 1 311 the 
contributions to the color factors 0(TpNfCA,F) of the Wilson coefficients L^ S ,H^ S ,H^ 
were also calculated. The corresponding Feynman diagrams consist of graphs with two 
internal fermion lines, out of which one is massless and one massive. After applying 
algebraic relations ||14|l these contributions to the massive Wilson coefficients can be 



represented in terms of the known set of weight w = 4 harmonic sums IT 1 311 



The O(aj) corrections were given in iH; see also the numerical implementation in Mellin space 0. 



Related cases concern the O(a^) contributions with two massive lines TpCp^ 
with either equal or different masses. In this note we present first results for these 
contributions. 



TWO MASSIVE QUARKS OF EQUAL MASS 

We first consider the case of two quarks of equal mass in the 0(a~) operator matrix el- 
ements oc Tp. In calculating the corresponding Feynman diagrams we apply the Mellin- 
Barnes representation [fl6ll through which the Feynman parameter integrals can be ex- 
pressed in terms of Meijer G-functions in general. In these calculations, like those 



in H3D, the use of modern summation algorithms as encoded in Sigma fl 1 811 are of es- 
sential importance. For the quarkonic flavor non-singlet case one obtains for the constant 
contribution to the operator matrix element : 
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The results in the flavor pure- singlet case are : 
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Here Pk(N) denote certain polynomials in N and = S$(N) are harmonic sums 111911 
These relations generalize the moments being obtained previously in [4]. The contribu 
tions to the anomalous dimensions given in P,!^] are confirmed. 



TWO MASSIVE QUARKS OF UNEQUAL MASS 

Beginning with O(ct^) graphs with internal massive fermion lines carrying unequal 
masses contribute. Since the mass ratio in case of the charm and bottom quarks is 
given by x = m%/ml ~ 1/10, one may expand the corresponding diagrams using this 
parameter. We first calculated the 2nd and 4th moment of the gluonic operator matrix 
element Ag g extending the code qexp [21] to higher moments applying projectors 



similar to those used in [4]. The unrenormalized operator matrix elements are given 
by: 
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where mi < m2 : m\ <C jj, 2 . The 4th moment reads : 
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These contributions to the massive Wilson coefficients can be uniquely calculated in 
the fixed flavor number scheme with three massless quarks in the initial state. However, 
they cannot be attributed either to the charm or bottom distribution in a variable flavor 
scheme. H This shows one of the limitations of this intention despite the fact that 
operator matrix elements are process independent quantities. In d3HU) power corrections 
contribute, however not of the kind 0(mj/Q 2 ) but of O^nf/m 2 ). The computation of 
the corresponding contributions for general values of N is in progress. 

We remark that the representation of the heavy flavor Wilson coefficients to O(a^) 
are given in the fixed flavor number scheme. As has been shown in Refs. 11231. 12411 this 
choice is sufficient for the kinematic region at HERA. 
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2 As has been shown in l22ll the application of a variable flavor scheme usually leaves no free choice of 
the matching scale jx 2 . In particular it is often quite different from the respective heavy flavor mass m?>. 
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